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Abstract 



Continuing our studies of radiative T decays, we report on a search for T — > 
77? and T 7/j(2220) in 61.3 pb"^ of e+e" data taken with the CLEO II 
detector at the Cornell Electron Storage Ring. For the 77/ search the three 
decays of the 77 meson to 7r~^iT~TT^, tt^tt^tt^, and 77 were investigated. We 
found no candidate events in the (Stt)*^ modes and no significant excess over 
expected backgrounds in the 77 mode to set a limit on the branching fraction 
of B{T 77/)<2.1 X 10"^ at 90% C.L. The three charged two-body final states 
hh {h = -K^,K'^,p) were investigated for /j(2220) production, with one, one, 
and two events found, respectively. Limits at 90% C.L. of I3(T — > 7/j) x 
B{fj hh) ~ 1.5 X 10~^ have been set for each of these modes. We compare 
our results to measurements of other radiative T decays, to measurements of 
radiative J/tp decays, and to theoretical predictions. 
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I. MOTIVATION 



Exclusive radiative decays of the heavy vector states J/ip and T have been the subject 
of many experimental and theoretical studies. For the experimenter, the final states from 
V 'jR are easy to identify and measure in that they have a high energy photon and 
low multiplicity of other particles. Backgrounds also tend to be small. Theoretically the 
emission of the photon leaves behind a glue-rich environment from which to learn about 
the formation of resonances from two gluons or to discover new forms of hadronic matter. 
Because gluons (g) themselves carry the quantum number of "color" , QCD allows for states 
with no valence quarks (g), but only gluon constituents: "glueballs". QCD also allows for 
more exotic combinations such as qgq "hybrids". These glueballs and hybrids are not just 
more resonances in the spectrum - they represent fundamentally different forms of matter 
from the mesons and baryons with which we are so familiar. 

For the J/ip charmonium system (the l^S^ state of cc) many such radiative two-body 
decays have been observed |Q, with some of the dominant ones being 7/2(1270), 777, and 
777'. However, the only such observation in the radiative decay of the T (the 138^ state of 
bb) is by CLEO in the final state 'yirn, in which an enhancement in the di-pion invariant 
mass consistent with being the /2(1270) meson was observed. A recent CLEO search ^ for 
the radiative production of the rj' meson yielded only an upper limit. 

In this Article we present a search for the radiative production of the other isoscalar 
pseudoscalar, namely the r] meson. As with the final state 777', this channel has received 
significant theoretical attention. The work of Korner, Kiihn, Krammer, and Schneider |^ and 
the followup publication by Kiihn ^ use highly virtual gluons to predict minimal suppression 
of radiative pseudoscalar production as the vector meson mass goes from that of the J/ ip to 
that of the T. Intemann used the vector dominance model (VDM) to predict branching 
fractions for T ^ 777 taking into account the interference between the T and T', the major 
contributing vector mesons in the model. Baier and Grozin showed that for lighter vector 
mesons (such as the J / ip) there might be an additional "anomaly" diagram that contributes 
significantly to the radiative decays. Ball, Frere, and Tytgat worked along similar lines [§. 
However, Baier and Grozin note that their approach applies directly to the "singlet" member 
of the meson nonet. Feldmann, Kroll, and Stech pursue the ideas of mixing in the decay 
constants of the pseudoscalars to derive ratios of their radiative production. Chao |]lOl has 
taken this approach further, determining mixing angles such as A^^^^ between the rjf, and rj in 



order to calculate radiative branching fractions. Finally, the recent work of Ma Jian-Ping |11 



uses factorization at tree level with non-relativistic QCD matrix elements to describe the 
heavy vector meson portion multiplied by a set of twist-2 and twist-3 gluonic distribution 
amplitudes. 

The other search we present here is for the radiative production of the /j(2220), also 
known as the .^(2230), in T decay. Many theoretical calculations of the spectrum of glueballs 
predict a J^*" = 2'^'^ state in the area of 2.2 GeV/c^. A candidate for this tensor glueball 
has been seen by some experiments, but not by others [^. The most complete claim of 
observation is by the BES collaboration who have published results [|12| for J/ip ^ 7/j(2220) 
with the /j(2220) reconstructed in tt+tt", K~^K~ , pp, and KgK^ as well as for 7r°7r° [|l3l and 
for 7777' ||T^. Of significant interest are several non-observations. CLEO has not seen the fj 



in two-photon interactions |15], which would lend credence to its being a glueball. However, 
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a narrow resonance in this mass region was not seen in pp production by either JETSET 
16| or, more recently, by Crystal Barrel [|l^]. This non-observation sheds doubt on the very 
existence of the fj. 



II. DETECTOR AND DATA SAMPLE 

Our analyses used 61.3 pb~^ of e~^e~ data recorded at the T(1S) resonance (^i = 9.46 
GeV) with the CLEO II detector [0 operating at the Cornell Electron Storage Ring (CESR). 
This corresponds to the production of A'^x = (1.45±0.03) x 10^ T(1S) mesons [§]. In addition, 
significantly larger samples taken near in time to this T(1S) data but at energies at or just 
below the T(4S) were used for comparison to the underlying continuum. The momenta and 
ionization loss [dE/dx] of charged tracks were measured in a six-layer straw-tube chamber, a 
ten-layer precision drift chamber, and a 51-layer main drift chamber, all operating in a 1.5 T 
solenoidal magnetic field. Photon detection and electron suppression were accomplished 
using the high-resolution electromagnetic calorimeter consisting of 7800 thallium-doped Csl 
crystals. The work presented here used only events with the primary, high-energy photon in 
the barrel portion of this detector, defined as | cos^| < 0.71, because the energy resolution 
for photons and reconstruction efficiency of the recoiling neutral mesons are degraded in 
the end cap regions and because the efficient, well understood triggers involve only the 
barrel region of the calorimeter. Between the central drift chamber and the electromagnetic 
calorimeter, strips of scintillating plastic were used for triggering and for measuring time- 
of-flight. Proportional tracking chambers for muon identification were located between and 
outside the iron slabs that provide the magnetic flux return. The Monte Carlo simulation of 



the detector response was based upon GEANT [|19[, and simulation events were processed 
in an identical fashion to data. 



III. SEARCH FOR T ^ 777 

Our search for T — > 777 involved the decays 77 77, 77 -k^-k^tx^, or 7/ — > tt^tt^tt'^; the 
latter two will collectively be referred to as (Bvr)". We followed procedures very similar to 
those used in our recent publication on the 777' flnal state In order to maximize detection 
efficiency and minimize possible systematic biases, we employed a minimal number of selec- 
tion criteria, with combinatoric background largely suppressed by requiring reconstruction 
of both the T and 7/ mesons. 

Events were required to have the proper number of quality tracks (either zero or two) of 
appropriate charges and at least three calorimeter energy clusters (which may or may not 
be associated with the tracks), of which one had to correspond to an energy of at least 4 
GeV and be in the barrel fiducial volume (| cos6'| < 0.71). In addition, we required that the 
events pass trigger criteria EO], based purely on the calorimeter, that were highly efficient 



and could be reliably simulated. 

For reconstructing vr'^ candidates, the photon candidates had to have minimum deposi- 
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tions of 30 (50) MeV in the barrel (endcap) regions^ and could not be associated with any- 
charged track; in addition, at least one of the two photons had to be in the barrel region. 
The 77 invariant mass had to be within 50 MeV/c^ (~ ±9o"7r) of the known 7r° mass [0; 
such candidates were then kinematically constrained to that mass. The photon candidates 
used in reconstructing the rj meson in 77 had to deposit a minimum of 60 (100) MeV in the 
barrel (endcap) calorimeter regions, could not be identified as a fragment of a charged track 
deposition, and had to have a lateral profile consistent with that of a photon. 

For the 7(87?)" modes we then built rj candidates from tt^'tt^tt'' or 7r"'"7r~7r°. Simulation 
events were used to determine the detector mass resolution for these two signal modes: 
o"^ = 10.7 and 9.0 MeV/c^, respectively. Candidates had to be within ±3cr^ of the known t] 
mass. In the case of the 77r°7r°7r° final state, no photon could be common to more than one 

combination. To suppress QED backgrounds in the 77r+7r^7r° final state, a charged track 
was rejected if its momentum, p, from the drift chamber matched its energy, E, as measured 
in the calorimeter as 0.85 <E/p< 1.05. 

Then, T candidates were formed by combining the high-energy photon (i?>4 GeV) with 
the r] candidate, requiring that this photon not already be used in reconstructing the event. 
To be considered, such a candidate had to have an invariant mass within ±300 MeV/c^ of 
y/s = rrir, a window of roughly three times the detector resolution as obtained from our 
simulations. Although, in general, multiple candidates per event were not restricted, there 
were two exceptions: (i) in the case of rj —>■ Tr^Tr^vr", if two rj candidates shared more than four 
photons, the candidate with the better combined for mass fits to the three vr" candidates 
was accepted; and (ii) in the case of rj —* vr^vr'vr'^, if two candidates for the neutral pion 
shared a daughter photon, the one with the better fit to the 7r° mass was taken. 

After these highly-efficient procedures were applied, we found no candidates in either the 
61.3 pb~^ of T(1S) data or in 189 pb~^ of continuum data samples.0 

From Monte Carlo simulations, the overall efficiencies for each channel, e^, were deter- 
mined to be (7.6 ± 0.8)% and (26.7 ± 1.5)% for the decay chains ending in rj tt^tt^tt^ 
and 1] —>■ TT^7c~TT^, respectively. The uncertainties here include the statistics of the Monte 
Carlo samples and our estimates on possible systematic biases, which we discuss be- 
low. Including the branching fractions for the rj decays |]1| and their uncertainties gave 
EhB^A = (8-7 ± 0.5)% for T ^ 7(37r)0. 

For the final state 777, the significant QED background compelled us to change the order 
of the constraints on the meson reconstruction and to add one additional selection criterion. 
Here we first took three photons, as defined above, and required that | I < 300 

MeV/c^. Then we took the photons in pairs and plotted the spectrum of m^^. 

A large background exists in both the T(1S) and T(4S) data sets, peaking near 0.40 
GeV/c^. From scanning such events it is evident that these are e'^e^ 77 events with pair 



The endcap region is defined as 0.85<| cos 0|<O.95; the region between this and the barrel fiducial 
region is not used due to its poor resolution. 

^ While some of these data were at the T(4S) resonance, we note that B meson decays cannot have 
a high energy photon so that the analysis of these data yields only udsc or QED backgrounds. 
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production by one of the photons giving a final state of 7e+e~. Tracks were not reconstructed 
in this subset of events due to the timing characteristics of the energy-based triggers in 
CLEO II. Such conversion events have their showers separated only in azimuth, in that the 
lepton pair has zero opening angle and any observed separation is due only to the bending 
by the magnetic field. Therefore, to suppress this 7e"'"e~ background, we removed events for 
which the angular separation in the calorimeter had no polar angle {6) component. 

To show that most of the remaining 777 background is of QED or continuum origin 
and not from the T(1S), we took the higher-energy data and performed the same analysis. 
We then subtracted this spectrum of m^^ from the T(1S) spectrum, after scaling it by 
the relative integrated luminosity, the relative reconstruction efficiencies, and the relative 
production rates. The last of these is determined from udsc continuum simulations at the 
two energies, for which we counted the number of events having one two-photon combination 
with an invariant mass within 3a^^ of m^. The mass resolution of cr^^ = 15.7 MeV/c^ is 
taken from simulation of signal events. The result of this subtraction is shown in Figure ^ 
the integral of the entries in this figure is 4.3 ± 9.4, consistent with zero. 
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FIG. 1. Shown as solid circles is the di-photon spectrum for T(1S) data after subtraction of 
the T(4S) data, appropriately scaled. The T(1S) data are fit, using a likelihood technique, for a 
T(1S)— > 77/ signal plus the scaled T(4S) spectrum, as detailed in the text. The superimposed curve 
shows the signal portion of the fit result. 



We then performed a binned likelihood fit of the T(1S) data, using the value of the 
corresponding bin in the scaled distribution from the T(4S) data and a Gaussian signal 
function. Here the Gaussian had a mean of the established value of the mass of the rj 
meson and a width taken from our simulations of the signal {a^j = 15.7 MeV/c^). The result 
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of the fit is shown in Figure |1|, yielding 4.0±3.8 events. Bin-by-bin statistical variations in the 
T(4S) spectrum were taken into account by performing this procedure multiple times, each 
time randomly assigning the number of events in each T(4S) bin, according to the statistics 
of that bin. The distribution function used in determining the limit was the average of these 
several functions. 

For this mode our overall reconstruction efficiency from signal simulation, including possi- 
ble systematic biases and statistical uncertainties from the simulation, is e^^ = (27.7±1.9)%. 

The major sources of possible systematic uncertainty in our efficiency calculation for 
T 'jr] are shown in Table |. These follow closely those detailed in our prior study of 
T 'J1]' 0. The degree of uniformity and definition of the fiducial volume of the barrel 
calorimeter (±2.2%) relates to our modeling the detector response to the proper angular 
distribution for this radiative T decay. Uncertainties in charged-track reconstruction (±1.0% 
per track), and trigger effects (±2.0%) were determined from previous detailed CLEO studies 
of low-multiplicity r-pair and 77 events. Similar studies allowed determination of the possible 
uncertainty in the reconstruction of 7r° and rj mesons from photons |2l| (±3% per meson); 
this relates to our ability to find and measure the daughter photons and, in the case of the 
neutral pions, to have the two-photon invariant mass be within the ±9(T7r window around the 
established 7r° mass. Our ability to model the E/p requirement in the 77r"'"7r~7r° final state 
was assessed using charged pions from Kg decays and assigned an uncertainty of ±3.2%. 
Shower leakage and other calorimeter effects make the mass distribution for T candidates 
asymmetric; based on CLEO experience with exclusive radiative B meson decays we 
have assigned an uncertainty of ±2% regarding our ability to model these effects. Based on 
a study of varying the mass resolution, a^^, in fits to the data before all the final criteria were 
imposed, we assign a systematic uncertainty of ±5% from this source. These uncertainties 
were added in quadrature, along with the statistical uncertainty associated with the size of 
Monte Carlo samples, to obtain the overall systematic uncertainty in the efficiencies. 

TABLE L Systematic uncertainty contributions, as relative percentages, to the efficiency for 
the studied decay modes for T 77/. Tlie combined uncertainties were obtained using quadrature 
addition. 



Uncertainty source 


77 




vr^vr vr*^ 


Fiducial requirements 


2.2 


2.2 


2.2 


Track reconstruction 






2.0 


7], TT^ reconstruction from 77 


3.0 


9.0 


3.0 


E/p criterion 






3.2 


Trigger simulation 


2.0 


2.0 


2.0 


T mass distribution 


2.0 


2.0 


2.0 


Variation of cj^^ in the fit 


5.0 






Monte Carlo statistics 


L6 


2.4 


1.0 


Combined uncertainty 


7.0 


10.0 


6.0 



The systematic uncertainties and the statistical uncertainty in the number of T(1S) 
decays are incorporated by a Monte Carlo procedure to obtain likelihood distributions for 
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the branching fraction in each mode as B{T —>■ '-yrj) = N^,/ (eA^x)- In this approach we produce 
multiple experiments with iV^ from the likelihood function appropriate for each decay modeQ 
and then divide by an efficiency and by a number of T(1S) events, each picked from a 
Gaussian distribution about their mean values with the appropriate standard deviations. 
Summing the resulting likelihood distributions, as shown in Figure 0, to 90% of their areas 
resulted in limits for 10^ x B{T ^ 777) of 28.2, 6.7, 2.6, 1.9, and 2.1 for 777, 77r°7r°7r°, 
77r"'"7r~7r°, 7(87?)°, and all combined, respectively. The number of observed events, detection 
efficiencies and limits are presented in Table 

TABLE II. Results for the search of T ^ 77/. Results include statistical and systematic uncer- 
tainties, as described in the text. 





77 






observed events 

Br^,^ (%) 

Reconstruction efficiency (%) 


4.0 ±3.8 
39.2 ±0.3 
27.7 ± 1.9 




32.2 ±0.4 
7.6 ±0.8 




23.1 ±0.5 
26.7 ± 1.5 


B{T ^ 77?) (90% C.L.) 


< 28.2 X 10-^ 


< 6.7 X 10^5 


< 2.6 X 10"^ 


Combined result 


< 2.1 X 10"^ 



To show that we could use CLEO data to observe 7r° and rj mesons we also applied our 
same selection criteria, with the exception of requiring a high energy photon, to samples 
taken at the T(1S) and at or near the T(4S). Figure ^ shows examples of inclusive yields of 
r] mesons in the T(4S) data that are consistent with the expected rates |l[]. 



IV. SEARCH FOR T ^ 7 Fj(2220) 



Our search for T — >■ 7/j(2220) used the three decay chains observed by BES ||T2| that 



involve two charged tracks: fj tt+tt^, fj K^K~ , and fj pp. We followed procedures 
very similar to those used in our publication on the observation of two-body radiative 
decays in T — 771 tt. 

Events were required to have two quality tracks of opposite charge with one energy 
deposition in excess of 4 GeV in the barrel fiducial volume (| cos6'| < 0.71) of the detector. 
The events were also required to pass trigger requirements pUf that were highly efficient for 
this process and that could be reliably simulated. 

Backgrounds from QED processes such as 7(7)yU"'"/i~ and 7(7)e"'"e~ are potentially large, 
so we next imposed criteria to minimize them. At least one of the charged tracks had to 
extrapolate to the barrel muon detector and have momentum above 1.0 GeV/c and yet not 
be identified as a muon. For each track the ratio of calorimeter energy to its measured 



^For the two 3(7r)'^ modes, which have no events actually observed, this likelihood function is a 
falling exponential as according to Poisson statistics. 



9 



1601001-014 




1 2 3 

Branching Fraction (10 ^) 

FIG. 2. Likelihood functions for branching fraction from the three final modes studied in our 
analysis and the combined likelihood function. All distributions include smearing by systematic 
uncertainties and have been normalized to unit area. The dotted vertical line is at 90% of the area 
of the combined function, namely 2.1 x 10^^. 



momentum had to satisfy E/p < 0.85 or E/p > 1.05, i.e., inconsistent with that expected 
for an electron. The tracks could not be consistent with coming from a photon conversion 
and had to have an opening angle between them of less than 162°. For each of the three 
decay modes [h = 7T'^,K'^,p) we established four-momentum conservation by demanding 
—0.03 < {E^fj^ — y/s) / < 0.02 and \p^j{fi\ < 150 MeV/c; these represent approximately 
±3(T in our detector resolution and take into account the effect of initial state radiation. 
We did not explicitly reject multiple combinations per event, and none were observed in 
our simulations of the signal. For each of the three modes we then plotted the di-hadron 
invariant mass, as shown in Figure ^. In that we will be comparing our results to those of 
BES [12|, we use a mass of rrif = 2.234 GeV/c^ and a natural width of Tf =17 MeV/c^ for 
the fj resonance. Our "signal box" is defined as ±2r/ centered at m/. For comparison, 
our signal simulations indicated mass resolutions (cr) of 15.6 ± 0.9, 14.8 ± 0.6, and 9.9 ± 0.4 
MeV/c^ for the final states vr^vr^, K^K~, and pp, respectively. 

The dominant continuum process ending in 77r+7r^ is e~^e~ 7p. We used this reaction in 
two ways. First, we repeated the prior analysis of two-body radiative decays to show that 
we could reproduce the shape and magnitude of the jp enhancement. To help us understand 
backgrounds we also studied data taken with an integrated luminosity of 838 pb~^ at energies 
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FIG. 3. The ^ 77 and rj tt^tt^tt^ invariant mass distributions from data taken at or near 
the T(4S). The plots give the invariant mass distributions (histograms), which are each fit with the 
sum (soHd hnes) of a polynomial background (dashed lines) and a Gaussian signal. 

at or just below the T(4S), recorded close in time to our T(1S) data. Before subtracting, 
we scaled the normalization of the tt+tt" invariant mass distribution from the higher energy 
data by C/{se) to take into account luminosity, the energy dependence of the cross section, 
and the reconstruction efficiency. We observe an excess of (29.4 ± 16.6) events for the T(1S) 
data in 0.45 GeV/c^ < mT,T, < 1.1 GeV/c^, consistent with zero excess. 

We studied backgrounds to the /j(2220) in three ways. When we determined upper 
limits we used the smallest of the three, thereby being the most conservative. In two of 
these, we used the data sets taken at or just below the T(4S) and accepted candidates 
within a background box width ilOFj centered at ruf. To understand the relative rates of 
production for background events at the two energies (^i = 9.46 vs. 10.56 GeV), we first 
simulated the udsc continuum, including initial state radiative effects, and for all events with 
a high-energy photon counted the number in which the invariant mass recoiling against the 
photon was within ±2r/ for the T(1S) energy, or within ztlOF/ for the higher energy;^ i.e., 
we did not fix the photon energy but rather the recoil mass region so that higher energy 
photons were involved in the events at i/i =10.56 GeV. The ratio of these numbers is 0.17 
as compared to the naive correction factor of 0.2 from the ratio of box sizes; i.e., the relative 
production rate at the T(1S) is 86% of that at the higher energy. 

If the dominant background in the signal region were from processes such as e'^e" 7p' 
or e+e~ 70, then the scaling would be as 1/s. In that case the T(1S) probability would 
be 124% that of the higher energy data. 

There may also be background contributions from other, as yet unmeasured, radiative 
T(1S) decays {e.g., T 7/4(2050)), which are not accounted for when comparing to the 
higher energy data. Therefore, as a third measure of the background, we look in the sidebands 



^This is similar to the approach used in the study of T ^ 777 — > 777 described in the preceding 
section. 



11 



of the fj region of the di-hadron spectra from the T(1S) data, namely 1.900 GeV/c^ < m^^ < 
2.200 GeV/c^ and 2.264 GeV/c^ < m,^^ < 2.500 GeV/c^. The number of events found in 
this region is then scaled by the ratio of the bin widths to predict the background in the 
signal region. 

For the channel fj tt^tt^, we observed one candidate in the ±2rj signal region of 
the T(1S) data set (see Figure |a). There were eight events in the broader ilOFj box at 
the higher energies, which when scaled by luminosity, efficiency, and the production ratio 
(as obtained from udsc simulations, described above) gave a mean background level in the 
signal box of fJ.bi'T^'^) = 0.12 events. Using the 1/s scaling or the sideband technique yielded 
0.17 and 0.36 events, respectively. 

Similarly for the channel fj —>■ K'^K~ we found one event in the signal box in the 
T(1S) data (see Figure ^). The analysis of the higher energy data showed 14 events in the 
iblOFj background region, distributed uniformly. When scaled appropriately for luminosity, 
efficiency, and production rate (as described above), this gave a mean background in the 
signal region of fib{KK) = 0.21 events. Using the 1/s scaling or the sideband technique 
yielded 0.30 and 0.72 events, respectively. 

The situation for fj pp was different. Assigning proton masses to the two charged 
tracks moves the large peak from e+e~ — >■ jp into the signal region. Although our selection 
criteria for conservation of four-momenta removed most of these backgrounds, we added a 
restriction that the time-of-flight as measured in the scintillation system be consistent for the 
two tracks to have proton masses, given their measured momenta. Also, nucleon-antinucleon 
annihilation can deposit large energies in the electromagnetic calorimeter, so our requirement 
on E/p would be very inefficient for anti-protons; therefore this requirement was removed 
for the negative track. The resulting distribution of pp invariant mass is shown in Figure ^jc, 
having two candidate events in the signal box region of ±2rj. In the higher energy data 
the background was no longer linear, due largely to the feed-through of e~^e~ •yp events. 
We therefore fit this distribution to a Gaussian for this 7p portion and a flat contribution. 
Integrating the fit result over the signal box region and correcting for relative efficiencies, 
integrated luminosities, and production (as described above) gave a mean background for 
this mode of Pbipv) = 0.28 events. A 1/s scaling would imply 0.40 events of background in 
the signal region, while the sideband technique gave a somewhat larger result of 0.48 events. 

In analyzing these three decay modes, we found that the signal regions are kinematically 
distinct so that no events are common to any two of them. 

To calculate reconstruction efficiencies we used our Monte Carlo simulations, based again 
on GEANT [|19|, with T — 7/j and fj hh for each of h = TT'^,K~^,p. The efficiencies. 



including possible systematic biases and uncertainties (as discussed below and summarized 
in Table 13) and the statistics of the Monte Carlo samples, were (28.8 ±4.2)%, (21.9 ±3.8)%, 
and (27.2 ±3.3)%, respectively. 

The simulation of signal events was generated uniformly in cos 6 of the high-energy pho- 
ton. Using the analysis for J/ip —>■ jfj with J = 2 ||2^ we have evaluated the range of 
possible angular distributions and their effect on our calculation of the geometric accep- 
tance, assigning a systematic uncertainty of ±10%. As in our search for T 777, we use 
previous detailed studies of low multiplicity t~^t~ and 77 events to study uncertainties in 
track reconstruction and momentum measurement (±1.0% per track) and trigger simula- 
tion (~ ±2.6%, and slightly dependent on momentum of the charged daughters). Also, as 
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before, we use previous CLEO experience with radiative B meson decays to assess pos- 



sible biases (±2%) in our T invariant mass determination due to shower leakage and other 
calorimeter effects. 

To evaluate the correctness of our simulations of the electron and muon rejection criteria, 
we analyzed charged pions from P^J^t with tt^tt^ (and charged conjugates) for 

both data and Monte Carlo. Using restrictive selection criteria on the masses of the 7r° and 
p~ mesons, on the specific ionization (dE/dx) of the charged pion, and the opening angle 
between the pions, we are left with a sample of charged tracks which is, according to the 
simulation, over 98% pions with less than 1% each of electrons and muons. Weighting the 
observed differences between simulation and data for these r data by the momentum spectra 
of the pions in the fj simulation shows that we need to increase our efficiencies by 2.6% and 
9.9% due to the effects of the E/p and muon requirements, respectively. We also assigned 
a systematic uncertainty of that same magnitude for these possible biases. For the final 
states K^K~ and pp, we scaled these uncertainties by the relative inefficiencies for events to 



pass these two criteria in our signal simulation, as indicated in Table |T|. We note that the 
magnitude of the uncertainty stemming from the muon rejection criteria is similar to that 
from our prior work [0 on T — > 77?"'" vr", which was ±13%, and calculated in an independent 
fashion. 

To estimate the systematic uncertainties associated with possible mismodeling of the 
time-of-flight identification we followed the strategy used in our study p4| of 77 pp. Here 



we varied the widths of the timing distributions by ±20%, which is about four times the 
precision to which they are known. This changes the reconstruction efficiency by ±2%, which 
we assigned as the systematic uncertainty from this source. 

TABLE III. Systematic uncertainty contributions, as relative percentages, to the efficiency for 
the studied decay modes for T — > 7/j(2220). The combined uncertainties were obtained using 
quadrature addition. 



Error source 


fj vr+vr" 


fj^ K+K- 


fj PP 


Angular distribution 


10.0% 


10.0% 


10.0% 


Trigger simulation 


2.6% 


2.6% 


2.7% 


Track reconstruction 


2.0% 


2.0% 


2.0% 


T mass distribution 


2.0% 


2.0% 


2.0% 


E/p criterion 


2.6% 


1.7% 


0.3% 


Muon suppression 


9.9% 


13.8% 


5.7% 


TOF identification 






2.0% 


Monte Carlo statistics 


1.0% 


1.2% 


0.9% 


Total (quadrature sum) 


14.8% 


17.6% 


12.3% 



To determine confidence limits for the results of this analysis, we followed the method 
advocated by Feldman and Cousins , which avoids under-coverage of confidence intervals. 
We adapted their method by replacing the confidence intervals for different mean numbers 
of observed events, fi, by confidence intervals for different values of the branching fraction, 
B = p/ (Nre). The uncertainties in the efficiencies and in the number of T(1S) produced were 
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incorporated by smearing the central values with Gaussian distributions. Finally, we extend 
the confidence limits so that they cover the integer values allowed by Poisson statistics. 

Before including systematic effects, we obtain 90% confidence intervals for the product 
B{T -ffj) X B{fj hh) of BB < 10.2 x 10"^ for tt+tt-, BB < 13.1 x 10"^ for K+R-, 
and the range 0.7 x 10^^ < BB < 14.3 x 10^^ for pp. Note that for the pp case we have an 
interval with both lower and upper limits. After adding these systematic effects, we obtain 
BB < 12.0 X 10-6 for tt+tt-, BB < 15.5 x 10-^ for K+R-, and the range 0.5 x 10"^ < BB < 
16.2 X 10"^ for pp. As noted earlier, we have used the estimate of the backgrounds, as 
determined by the simulated production ratios. Using larger background estimates, such as 
from 1/ s scaling or sideband comparisons, would give lower upper limits and eliminate the 
lower limit in the case of fj — > pp\ for these three background estimates the probability of 
observing two or more events is between 3 and 8%, insufficiently small to claim a signal in 
this decay mode. We therefore take the conservative approach and quote the 90% confidence 
level upper limits as given and quote no lower limit for any of the channels. The results of 
our analysis are summarized in Table 

TABLE IV. Results of the search for T — > 7/j. The estimated background, obtained by scahng 
by the simulated continuum production rates, is the smallest of the three estimates of background 
levels and leads to the most conservative upper limits. The efficiencies and limits include systematic 
effects. The last row of entries is a scaling of the BES results for J /tp decays. 





fj TF+TT- 


fj^ K+K- 


fj PP 


Observed events in ±2rj 


1 


1 


2 


Scaled continuum background (^5) 


0.12 


0.21 


0.28 


Overall efficiency 


(28.8 ±4.2)% 


(21.9 ±3.8)% 


(27.2 ± 3.3)% 










e(/j(2220) ^ hh) (90% C.L.) 


< 12.0 X 10-6 


< 15.5 X 10-6 


< 16.2 X 10-6 


B{fj{2220) hh) x 0.04 


1.5l0| X 10-6 


2.bt\'i X 10"^ 


0.7 ±0.3 X 10-6 



As a comparison we have also computed the upper limits following the traditional ap- 
proach of summing up the Poisson distribution V{< n\ fj,) = e~^x (1+/^+... + ^) to 90% C.L. 
For n = 1 and 2, this gives limits of /x = 3.88 and 5.32, respectively. Ignoring backgrounds 
the corresponding upper limits on the product branching fractions, including systematic ef- 
fects, are BB < 9.3 x 10-6 for tt+tt", BB < 12.2 x 10-6 for K+R-, and BB < 13.5 x 10-6 for 
pp. Note that this method can lead to under-coverage, as outlined in the paper by Feldman 
and Cousins 



V. COMPARISON TO PRIOR RESULTS AND THEORY 



The only other reported analysis [|T| of B{T ^rf) is by Crystal Ball [^, which deter- 
mined a 90% C.L. upper limit of 3.5 x lO-'^; our limit of 2.1 x 10"^ is ~ 17 times more 
stringent. 
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FIG. 4. Di-hadron invariant mass distributions for data events passing all selection criteria. 
The inserts show magnifications of the signal region. The vertical arrows in the inserts indicate the 
limits of the signal box, which are at ib2r, F = 17 MeV/c^. No events are common to the signal 
regions of these channels. 



We can then use our new limit on B{T 'yrj), the measured enhancement [0 near 1270 
MeV/c^ in T ^ 77r7r, and the measurements of BiJ/^p 7/2(1270))= (1.38 ± 0.14) x lO^^ 
and B{J/ip — > •yrj) = (0.86 ± 0.08) x 10~^ |Q to create interesting ratios. Here we assume 
that the observed structure in the di-pion mass spectrum near 1270 MeV/c^ in the T study 
is totally due to production, which implies B{T 7/2(1270))= (8.2 ± 3.6) x 10"^ We 
then obtain 



Rf{T,J/ij) 



g(T ^ 7/2) 
B(j/^ ^ 7/2 



0.061 ± 0.026 and 



" B^i?/^~^-Jt) ^ 

and, looking instead at the ratios of the two final states for a given vector parent. 



(2) 
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Rriv, /2) ^ ^ ^"I^, < 0.32 (90% C.L.) and (3) 
^( i ^ 7/2) 

R.,AvJ.)^ '^lf/^::;jl -o.e2±om. (4) 

These results were obtained by including the various uncertainties in a Monte Carlo tech- 
nique. In the last two of the ratios, we assume that all the uncertainties from the J/ip and 
T measurements are uncorrelated. Our limits show that the branching fractions into 777 and 
7/2(1270) behave differently in the cases of J/?/' and T, although not as dramatically as in 
the case of T — > 777' 0. In addition we form the double ratio: 

Rr{v,h) _ BjT^^v) g(J/^->7/2) , . 

^ - Rj/4v, h) ~ B{T 7/2) B{J/^ - ^v) < "-^^ 

at 90% C.L. This is to be compared with the prediction of Korner et al. 0] of 

7^^w, = ^ = 0-42; (6) 
Chao's technique [jlO| first calculates mixing angles among the various pseudoscalars, 



extending the J^^ = 0~"*" nonet to include heavier cousins such as the rjb. Then, using the 
predicted allowed Ml transition T — 777^, Chao predicts B{T — »• 777) = 1 x 10~^, which is 
consistent with our limit. We note that in our prior work p] we did not know of Chao's 
prediction of B{T — > 777') = 6 x 10"^, which is to be compared with our upper limit for that 
process of 1.6 x 10"^ at 90% C.L. 

Intemann's extended vector dominance model gives 6.5 x 10~^ < B{T '-frj) < 1.2 x 
10^'', with the two limits determined by having destructive or constructive interference, 
respectively, between the terms involving T and T'. This range of predictions is well below 
our new limit. Ma uses a technique in which the decay amplitude factorizes into a non- 
relativistic piece describing the bound state of the heavy quarkonium and an expansion in 
"twist" to characterize the conversion of the gluons into the final state meson. The published 
value [0 was B(T — > ■yrj) = 1.2 x 10~^, although subsequent correspondence |^ indicates 



the correct value is actually 3.3 x 10^^. In either case, these are significantly below our limit. 
Feldmann et al. predict the ratio r(T '-fi]') /r{T '~fi]) = 6.5. Given that we only 
have limits on these two processes, we cannot address their prediction. 

For the glueball candidate, our 90% C.L. limits on the product branching fractions B{T 
7/j) X B{fj hh) are on the order of 1.5 x 10~^ for each of the three modes. In Table |I^ we 



show the results from BES ||12| for these channels in radiative J/ ip decay, scaled by a factor 
of JF. This scaling arises from the naive expectation that the amplitude for the radiative 
process of meson formation varies directly as the quark charge (to couple to the photon) and 
inversely as the quark mass (from the fermionic propagator between the photon emission 
and the resonance formation). One then squares this to get the rate and corrects for the full 
widths of the heavy quarkonia to obtain 

^=(^)2.^^ = 0.04. (7) 
qcrrib Tr 
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Here we have used masses of 1.7 GeV/c^ and 5.2 GeV/c^ for the charm and bottom 
quarks, respectively. As shown in Table |I^ and in Figure ^, our limits on radiative /j(2220) 
production do not confront these predictions. 

As a quantification of this situation, we have have used a Monte Carlo technique to 
evaluate the ratio of the CLEO and BES |]T2| results, i.e., 



e(J/'* ^T/,;) B{J/t ^ lij) ■ B(fj ^ hh) 

For the denominator we add the BES statistical and systematic uncertainties in quadra- 
ture and throw BES "experiments" in a Gaussian distribution, keeping only physical values. 
In all three cases {h = vr"*", K~^,p) we can only say that Rfj is less than roughly 0.50 at 90% 
C.L., an order of magnitude from the predicted ratio of ^ = 0.04. 

We also note that the Crystal Barrel collaboration has combined their results with 



those of BES [H to obtain that BiJ/-^ 7/j(2220)) > 0.003 at 95% C.L. Applying the 
naive scaling factor would then predict that B{T 7/j(2220)) > 1.2 x 10""^; this would 
be larger than the CLEO result for the radiative decay to /2(1270). 

VI. SUMMARY AND ACKNOWLEDGMENTS 

In summary, we have used the CLEO detector operating at the CESR storage ring to 
search for two-body radiative T(1S) decays. In the work presented in this Article, we reported 
specifically on searches for the decay T — >■ 777 with the subsequent decays 77 tt'^tt^tt^, 
rj — >• 7r"'"7r~7r°, and 77 — >• 77 and for the decay T 7/j(2220) with the subsequent decays 
of the glueball candidate of fj tt+tt^, fj — > K^K^, and fj pp. Including our prior 
pubhshed results we have the following: 

for the decay T •ymr [||, 

B{T 77r+7r") = (6.3 ± 1.2 ± 1.3) x 10"^ [m^^ > 1.0 GeV/c^] , 

and 

B{T 77r°7r°) = (1.7 ± 0.6 ± 0.3) x 10~^ [m^^ > 1.0 GeV/c^] ; 

and, assuming the enhancement in the invariant mass spectrum for tt+tt^ in the region of 
1.3 GeV/c^ is due to /2(1270) production, 

B{T 7/2(1270)) = (8.2 ± 3.6) X 10'^; 
for the decay T ■jrj' 0, 

B{T -fT]') < 1.6 X 10"^ [90%C.L.]; 

for the decay T ^ 777, 

B{T 77/) < 2.1 X 10"^ [90%C.L.]; 
and for the product branching fractions involving the glueball candidate /j(2220), 

B{T -> 7/j) X B{fj 71+71-) < 1.2 X 10"^ [90%C.L.], 
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B{T ^ 7/j) X B{fj K+K-) < 1.6 x 10" 



C.L.l 



and 



C.L. 



B{T ^ 7/,,) X B{fj ^ pp)< 1.6 X 10-^ 

We also present this summary in graphical form in Figure |^: (i) the world-average [|l[ 
values for the radiative two-body decays of the J/ip to t], rj', and /2(1270); {ii) the BES 
12 1 results for the product branching fractions B{J/iIj — » ifj) x B{fj 



h+h- 



[iii] the 



scaling of these J/ip results by to obtain estimates for the corresponding T(1S) decays; 
{iv) the CLEO results summarized above; and [v) theoretical predictions from Korner et al. 
0], Intemann 0], Ma [ill], and Chao p[| . 
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FIG. 5. Radiative decays of J and T(1S) vector mesons. Shown are the experimental results 
(as solid circles or limit bars) from the PDG for the well-established radiative J/ip decays to ry, r//, 
and /2(1270), from BES for the three charged modes of radiative decay to the glueball candidate 
/j(2220), and from CLEO for the T decays. The solid triangles give the values for radiative T 
decay based on radiative J /ip decay and the naive scaling involving the masses and charges of the 
constituent quarks. Various explicit theoretical predictions for the radiative production of ?7, 77', 
and /2(1270) are shown as open symbols. For references see the text. 
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